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Abstract 

We present a detailed analysis of the Landau-Zener problem for 
an interacting Bose-Einstein condensate in a time-varying double- 
well trap, especially focussing on the relation between the full many- 
particle problem and the mean-field approximation. Due to the non- 
linear self-interaction a dynamical instability occurs, which leads to 
a breakdown of adiabaticity condition and thus fundamentally alters 
the dynamics. It is shown that essentially all features of the Landau- 
Zener problem including the depletion of the condensate mode can be 
already understood within a semiclassical phase space picture. In par- 
ticular, this treatment resolves the formerly imputed incommut ability 
of the adiabatic and semiclassical limits. The possibility to exploit 
Landau-Zener sweeps to generate squeezed states for spectroscopic 
tasks is analysed in detail. Moreover, we study the influence of phase 
noise and propose a Landau-Zener sweep as a sensitive, yet readily im- 
plementable probe for decoherence, since this has a significant effect 
on the transition rate for slow parameter variations. 
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1 Introduction 



The Landau-Zener problem aims at the general description of nonadiabatic 
transitions at avoided level crossings. In the standard setting, the dynamics 
is restricted to two levels with a constant coupling J, whose energy difference 
varies linearly in time, e(t) = at. Of particular interest is the Landau-Zener 
tunneling probability between the two adiabatic states which is found to be 

Plz = e- j2 / Q , (1) 

independent of the initially occupied level. Due to its generality, this result 
has been applied to numerous problems in various contexts like, e.g., spin- 
flip processes in nano-scale systems [T], molecular collisions [2], quantum- 
dot arrays [3], dissipative systems [HE] or quantum information processing 
tasks |&H9] , to name but a few examples. 

The Landau-Zener szenario was one of the first major problems adressed 
within time- dependant quantum theory. While the single-particle case was 
solved independently by Landau, Zener, Majorana and Stiickelberg already 
in 1932 [TUHT3] . the generalization of the results to interacting many-particle 
systems remains an open question up to today, and even the mean-field 
dynamics is not yet fully understood. The non-linear self-interaction fun- 
damentally alters the dynamics, leading to a breakdown of adiabaticity due 
to the bifurcation of nonlinear stationary states [T1H2T)] . The many-particle 
Landau-Zener problem is of fundamental interest not only from the theoret- 
ical but also from the experimental point of view and has in recent years 
attracted a lot of interest, especially in the context of the dynamics of Bose- 
Einstein condensates (BECs) in optical lattices [2TM25] . 

In the following, we present a detailed analysis of nonlinear Landau-Zener 
tunneling between two modes focussing on the relation between the original 
full many-particle problem and the mean-field approximation. The break- 
down of adiabaticity is a consequence of a bifurcation of the mean-field sta- 
tionary states or the occurrence of near-degenerate avoided crossings in the 
many-particle spectrum, which are intimately related. Furthermore we dis- 
cuss the Landau-Zener problem within a semiclassical phase space picture, 
where the quantum dynamics is approximated by a Liouvillian flow rather 
than a single trajectory. It is shown that essentially all features of the dy- 
namics including the depletion of the condensate mode can be already un- 
derstood within this approach. In particular, this treatment resolves the 
formerly imputed incommutability of the adiabatic and semiclassical limits. 



2 



Number squeezing effects during the transition are analysed in detail. More- 
over, we study the influence of phase noise which is an unavoidable feature 
in every experiment. Since it has a significant effect on the transition rate 
for slow parameter variations, a Landau-Zener sweep is a sensitive probe for 
decoherence. 

In particular, this paper is organized as follows: In section [2] we first intro- 
duce the many-particle and mean-field description of the system and define 
the Landau-Zener transition probability Plz in both cases. The basic fea- 
tures of the nonlinear Landau-Zener problem are reviewed in section [3j In the 
mean-field approximation, Plz does not vanish even in the limit a — > if the 
interaction strength exceeds a critical value. In contrast, the many-particle 
Landau-Zener tunneling probability always tends to zero in the adiabatic 
limit. However, this convergence is extremely slow so that the breakdown of 
adiabaticity is approximately present also in the many-particle description. 
To highlight the origin of this breakdown we analyze the full quantum state 
during a Landau-Zener sweep in more detail in section H] using the SU(2)- 
phase space techniques derived in [2S1I2Z]- We show that many features of 
the many-particle dynamics can be captured to astonishing accuracy within 
the phase space description, including the depletion of the condensate mode 
as well as number squeezing of the final state. Yet, we show that using a 
Landau-Zener sweep to generate squeezed states for quantum metrology is 
very difficult for realistic systems. Section then gives a detailed analysis of 
the region of validity of the mean-field approximation and the convergence to 
the mean-field limit. In section [6] we will briefly discuss the influence of phase 
noise, which is unavoidable in every experimental realization. We conclude 
with a short summary and outlook. 

2 Mean-field and many- particle description 
of a two-mode BEC 

The Bose-Hubbard type hamiltonian 

H = e(i) (4^2 - a\a^j - J (a\a 2 + a\a^j + ^ {h-i(h\ - 1) + h 2 {h 2 - 1)) (2) 

describes the dynamics of ultracold atoms in a double-well potential or the 
dynamics of a system of two-level atoms, respectively [28H3Tj . The operators 
hi and a 2 annihilate an atom in the first and second mode, respectively, while 
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the operators hi = a\a,i describe the population of the wells % = 1,2. The 
tunneling matrix element and the on-site interaction strength are denoted 
by J and U and the time-dependent energy offset of the two modes is given 
by 2e(t) = 2at. In all numerical examples we shall set H — J — 1, thus 
measuring time in units of the tunneling time h/J. 

The time evolution generated by the hamiltonian 02J preserves the total 
particle number. It can be rewritten using the generalized angular momen- 
tum operators 



L x = i (a\a 2 + 4«i 



Ly = o (al«2 - > (3) 



0^0*2 — Ct-]CL\ 



which correspond to the tunneling, the momentum and the population im- 
balance, respectively. In this this representation, the hamiltonian ([2]) is given 
by 

H = 2e(t)L z - 2JL X + UL 2 Z . (4) 

Initially, the two modes are energetically well separated and the ground state 
of the Bose-Hubbard hamiltonian (T5]) is 

|vI/(t^-oo)) = (iV!)- 1 /2(at 1 )^|0), (5) 

thus we assume that initially all particles are localized in the first well, corre- 
sponding to a fully condensed state. The many-particle Landau-Zener tran- 
sition probability for the population is then given by 

_ (nig -> +oo)) 

^ - (ril(t ^_ oo)) - w 

In the following, the many-particle quantum state is denoted by a capital 
while the lower case ip is used for the components of the mean-field state 
vector. To distinguish the transition probabilities we use the superscripts mp 
and mf for the many-particle and mean-field quantities, respectively. 

In the mean-field approximaton, the time evolution is given by the discrete 
Gross-Pitaevskii equation [30ll3Tj . 

il f ^ \ = ( at + 9\^i\ 2 -J \(^\ ( 7 ) 

dtlVW 1 - J -at + g\^ 2 \ 2 J\ ^ J ' {<) 
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where g = UN is the macroscopic interaction strength. The mean-field 
approximation is valid in the limit iV — > oo, while g is kept constant. In 
close analogy to the angular momentum operators fl3]), we define the Bloch 
vector 

S * = \ (^1^2 + ^2^1) , 

^ = ^(^2-^1). (8) 

The dynamics of the Bloch vector s is restricted to the surface of the Bloch 
sphere, as the norm s 2 = 1/2 is conserved by the equations of motion in 
the absence of phase noise. Thus, a convenient representation of the Bloch 
vector is given by the polar decomposition 



sin 9 cos 1 
9 sin < 
cos 9 



- J sin 9 sin ] . ( ) 



In this setting, the Landau-Zener tunneling probability in the level j = 1,2 
is defined as 

LZ '"hM^-oc)r (10) 

Again we assume that all particles are initially localized in one of the modes, 
i.e. ipj{t — > —00) = 1. 

A significant extension of the applicability of the mean-field approxima- 
tion is achieved if one considers the dynamics of quantum phase space distri- 
butions instead of single mean-field trajectories. While the common mean- 
field approach allows only statements about expectation values, the phase 
space description takes also the higher moments and their time evolution 
approximately into account. Here, we will only review the basic definitions. 
For further details and a rigourous mathematical introduction see [26|[27] and 
references therein. 

The starting point is the Husimi or Q-Function, which is definded as the 
projection onto the set of SU (2)-coherent states 

Q(M,t) = |(M|v|/(t)>| 2 , (li) 
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with 



1 / i\ N 
\6, 4>) = —= [cos{9/2)a\ + sin(0/2)e-^4J |0, 0). (12) 



_'N\ 

Note that the quantum state can be uniquely reconstructed from this repre- 
sentation, due to the overcompleteness of this basis set. The exact dynamics 
of the Husimi function is then given by 



| W, *) = {&(*) A + 2 J ( sin ^ - cos cot » £ 

-Qcose— + ^ S me 1 ^}Q(t).o). (13) 



It is important to note that this exact evolution equations can be written 
as a classical Liouvillian phase space flow plus a quantum correction term 
which vanishes as 1/N. The classical part is equivalent to the discrete Gross- 
Pitaevskii equation (j7J) in the appropriate parametrization [27]. A semiclas- 
sical approximation of the phase space flow is thus provided by a truncated 
phase space dynamics: The initial state is mapped to its Husimi function, 
which is then propagated according to a classical Liouville equation omitting 
the quantum corrections in equation ( 1T31) . Equivalently we will consider an 
ensemble of classical phase space trajectories whose starting points are dis- 
tributed according to the initial Husimi function. The truncated phase space 
evolution defined above clearly goes beyond the common mean-field dynam- 
ics as it enables us to approximate the dynamics of the higher moments of 
the quantum state. 

Within the phase space description, the expectation values of the gener- 
alized angular momentum operators (J3j) are obtained by an integration over 
the quasi-probability density: 

(L) = (N + 2) J s(0, <p)Q(9, 0) sin 9d6d(j). (14) 

As a direct consequence, we can calculate the reduced single particle density 
matrix (SPDM), which is defined as 

. = ( 1/2 - (L,)/N (L x )/N - i(L y )/N \ 

P \(L x )/N + i(L y )/N 1/2 + (L x ) /N J' 1 1 
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Figure 1: Landau-Zener tunneling probability in the lower level (left) and the 
upper level (right) for J = 1 and g = 0,1, 5. Mean-field results P^(oc) are 
plotted as solid blue lines, many-particle results P^P(a) for N = 50 particles 
as red circles. 

The SPDM is a very useful quantity, since it characterizes the many-body 
quantum state of the trapped atoms. In particular, the fraction of atoms 
condensed to a single quantum state (the BEC) is given by the leading eigen- 
value of the SPDM |32j. If the the expectation value (L) /N lies on the Bloch 
sphere, i.e. has a magnitude of 1/2 (as it is always the case in the com- 
mon single-trajectory mean-field approach), then the two eigenvalues of the 
SPDM are always {0, 1} indicating a pure BEC. The phase space represen- 
tation is not limited to product states. Due to the averaging procedure in 
equation f JT4|) . the expectation value of the Bloch vector is then no longer 
restricted to the surface, but can lie anywhere inside the Bloch sphere. The 
phase space approach has been proven to be a very useful tool to go beyond 
the usual mean-field description [27] , especially in the description of dynam- 
ical instabilities, where it is clearly not sufficient to take into account only 
expectation values and to neglect all higher moments. As we show in the 
following, this approach also resolves the non-commutability of the adiabatic 
and the semiclassical limit, which therefore must be considered as an artifact 
of the single-trajectory description. 

3 Nonlinear Landau-Zener tunneling 

The nonlinear self-interaction fundamentally alters the dynamics of the sys- 
tem [Tll[T5"| IT7] and strongly influences the Landau-Zener transition probabil- 
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Figure 2: Energy of the mean- field stationary states ( fl6i) in comparison to 
the eigenenergies of the two-mode Bose-Hubbard model ([2]) in dependence of 
the offset e for J = 1 and g = 1 (left) and g = 5 (right) and iV = 20 particles. 

ity, as can be seen in figure [TJ The solid lines show the mean-field Landau- 
Zener tunneling probability ( flOl) in dependence of the parameter velocity 
a for different values of the interaction strength g. For this calculation we 
have used the common single-trajectory mean-field approximation. However, 
there are no visible differences to the phase-space results for the actual pa- 
rameters. The open circles represent the corresponding many-particle results. 
In the linear case g = 0, one recovers the result ([1]) for the Landau-Zener 
tunneling probability. For a slow parameter variation, the state can adia- 
batically follow the instantaneous eigenstates and thus most particles tunnel 
coherently to the other well. For a faster sweep, this coherent tunneling effect 
is strongly disturbed such that the Landau-Zener transition probability no 
longer vanishes. This effect is present in both, the transition probability in 
the upper and the lower level. 

In the nonlinear case, the tunneling probability becomes strongly asym- 
metric: it increases as g increases in the upper level, while it decreases in 
the lower level. To understand this effect, it is insightful to consider the 
total energy of the mean-field system. Figure [2] shows the eigenenergies of 
the Hamiltonian ([2]) in comparison to the total energies of the 'nonlinear 
eigenstates', i.e. the stationary states of the mean-field dynamics (J7J), 

E°* = e(|^ 2 | 2 - |^| 2 ) - J(^*V>2 + tfttk) + f (l^il 4 + |^ 2 | 4 ). (16) 

Compared to the non-interacting case g = 0, the left-hand side shows that 
the upper level is sharpended, while the lower level is flatened for small 
interactions < g < 2J. This flattening supresses the tunneling probability 
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Figure 3: Dynamics of the many-particle Landau-Zener transition for J = 1, 
g = 5, iV = 20 and a = 0.1. The upper figures show a contour plot of 
the Husimi distribution of the instantaneous eigenstates marked in the level 
scheme on the right at times tj = —10,0.5,5.5,16. The lower figures show 
the Husimi distribution of the dynamical quantum state at the same 

times. 

from the lower level to the upper level, leading to a decreased Landau-Zener 
probability in the adiabatic regime. On the other hand, the sharpening of 
the upper level makes it more difficult to follow the adiabatic eigenstates, 
which results in a increased Landau-Zener probability for the upper level, as 
can be seen on the right-hand side of figure [TJ 

Most remarkably, the tunneling probability in the upper level does not 
even vanish in the adiabatic limit a — > for g > 2J, i.e. adiabaticity breaks 
down in the strongly interacting case. 

In order to explore the origin of this breakdown of adiabaticity we com- 
pare again the eigenstates of the many-particle system to the stationary 
states of the mean-field system. For g > 2J the mean-field eigenenergies 
show a swallow's tail structure in the upper level, reflecting the occurrence 
of a bifurcation of one of the steady states into three new ones, one of them 
hyperbolically unstable (dashed line) and two elliptically stable (solid lines). 
The system can adiabatically follow the steady states as long as these are el- 
liptically stable. This is possible only until the end of the swallow's tail where 
the elliptic fixed point vanishes in an inverse bifurcation with the hyperbolic 
fixed point [TSJEIS]- Then the dynamics becomes unstable and adiabaticity 
is lost even for very small values of a. 
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The swallow's tail in the mean-field energy corresponds to a caustic of the 
many-particle eigenenergy curves in the limit N — > oo, which are bounded by 
the mean-field energies from below and above. Within this caustic one finds a 
series of quasi-degenerate avoided crossings of the many-particle levels. The 
level splitting at these crossings tends to zero exponentially fast in the mean- 
field limit N — > oo with g = UN fixed [TSJEl] . Thus the system will show a 
complete diabatic time evolution at these quasi-crossings even for very small 
values of a. Outside the swallow's tail one finds common avoided crossings, 
where the system evolves adiabatically for small value of a. 

Note, however, that the breakdown of adiabaticity is only approximate 
for the many-particle system. It is known that for a symmetric tridiagonal 
Hamiltonian, such as the one we are considering (T5]) with J ^ 0, the level 
spacings in the spectrum may be exponentially small but nevertheless always 
non-zero [33] . Thus adiabaticity can be restored when the parameter velocity 
a is decreased well below the square of the residual level splitting A [34|: 

i 

a < A 2 with A oc Nexp(-r)N), (17) 

where rj is a proportionality constant which depends of the system parame- 
ters. However, the adiabaticity condition on the velocity becomes exponen- 
tially difficult to fulfill. Thus the breakdown of adiabaticity is also present in 
the full many-particle system for any realistic set of parameters. The same 
dynamics is found for attractive nonlinearities, g < 0, only the roles of the 
upper and lower level are exchanged. 

4 Phase space picture 

Further insight into the dynamics of nonlinear Landau-Zener tunneling can be 
gained within the phase space picture introduced in section [2] and in ffiH ITr] . 

According to the remarks in the previous section, the system will un- 
dergo a series of diabatic transitions up to the end of the swallow's tail and 
evolve adiabatically afterwards. To verify these claims, we compare the ac- 
tual many-particle quantum state to the instantaneous eigenstate in 
figure [3] at four points in time during a Landau-Zener passage. To visualize 
the quantum states, we use the Husimi distribution Q(9,(f>,t) as defined in 
equation ffTTj) . The right-hand side of figure [3] illustrates the series of dia- 
batic/adiabatic transitions and the specific instantaneous eigenstates shown 
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Figure 4: Dynamics of the many-particle Landau-Zener transition for g = 5, 
N = 50 and a = 0.1. The left figure shows the eigenvalues of the SPDM 
( fl~5|) . The expectation values (£/&) and widths AL^ of the angular momentum 
operators ^ for k — x (blue), k = y (red) and k = z (green) are plotted 
on the right-hand side. The onset of the oscillatory dynamics corresponds 
to the cusp of the swallow-tail structure in the mean-field representation (cf. 
figure El) . 

in upper panels of the figure. One observes a good agreement between the 
dynamical state and the instantaneous eigenstates, during the transition as 
well as afterwards. However, the crossover from diabatic to adiabatic tran- 
sitions is not absolutely sharp. The final state contains small contributions 
from other instantaneous eigenstates. 

In order to characterize the many-particle quantum state during the 
Landau-Zener transition, we have plotted the eigenvalues of the SPDM ( 1X5]) 
on the left-hand side of figure HI One eigenvalue remains equal to unity, while 
the other one vanishes, indicating a fully coherent state until the crossover 
from diabatic to adiabatic transitions. Then one observes an oscillation of the 
SPDM eigenvalues: The contributions of the different many-particle eigen- 
states de- and rephase periodically giving rise to a beat signal which is gen- 
uinely quantum. The oscillation of the coherence is mirrored in the evolution 
of the uncertainties of the angular momentum operators AL X and AL y shown 
on the right-hand side of figure HI The uncertainties are strongly enhanced 
when the coherence is (partly) lost. This behaviour can be intuitively ex- 
plained in terms of the dynamics of the Husimi distribution. The centre of 
mass of the Husimi function oscillates rapidly in the 0-direction, leading to 
oscillations of the expectation values (L x ) and (L y ). Furthermore the dis- 
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Figure 5: Semiclassical simulation of the many-particle Landau-Zener tran- 
sition for J — 1, g — 5, N = 50 and a = 0.1. The left figure shows the 
eigenvalues of the ensemble approximation for the SPDM = (ipl^i) ci- 
The ensemble expectation values N(sk) c \ an d the standard deviation NAs^ 
of the Bloch vector (jSJ) for k = x (blue), k = y (red) and k = z (green) are 
plotted on the right-hand side. 

tribution breathes in the ^-direction at a slower timescale, leading to the 
oscillations of the width AL X and AL y and the periodic revivals of the co- 
herence. The oscillations of the expectation values die out at the times when 
the Husimi function is spread nearly uniformly in the ^-direction, i.e. at 
the times where the coherence is minimal. In contrast, the Husimi distribu- 
tion is well localized in the ^-direction for long times and the corresponding 
uncertainty AL Z remains small. The population difference (L z ) is thus well 
described by the simple Bogoliubov mean-field approximation. Many-particle 
and mean-field results for the Landau-Zener tunneling rate show an excel- 
lent agreement (cf. figured]), because they depend only on the population 
difference and not on the coherence. 

The evolution of the coherence and the uncertainties AL X and AL y cer- 
tainly goes beyond the Bogoliubov mean-field approximation, but most of the 
effects can be taken into account by the semi-classical phase space approach 
introduced in section [2] and in [26j|27]. Figure |6] shows the dynamics of the 
many-particle Landau-Zener scenario in quantum phase space in comparison 
to the dynamics of a classical phase space ensemble. The expectation values 
and variances of the Bloch vector L calculated from such an ensemble simu- 
lation are plotted in figure |5j It is observed, that the spreading of the Husimi 
distribution in the direction of the relative phase 4> and the loss of coherence 
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Figure 6: The many-particle Landau-Zener scenario in phase space. The 
dynamics of an ensemble of 150 classical trajectories (upper panels) is com- 
pared to the evolution of the Husimi distribution (lower panels) at times 
t — —16, —8, 0, 8, 16 (from left to right). Parameters are chosen as in figure 

El 

are well reproduced by the classical ensemble. However, the quantum beat 
oscillations of the coherence are of course not present in the classical dis- 
tributions as shown in figure [5j The expectation value and the fluctuations 
of the classical Bloch vector s defined in equation ([8]) show a similar effect. 
The global dynamics of the angular momentum operator L plotted in figure 
H]is well reproduced, whereas all the quantum beats are absent. These are 
genuine many-particle quantum effects. 

The previous results show that the many-particle quantum state after a 
nonlinear Landau-Zener sweep is far from being a pure BEC. In particular 
it has been claimed that the final state is strongly number squeezed in com- 
parison to a pure BEC with the same density distribution [12]. The figures H] 
and show the evolution of the expectation values and variances of L, com- 
paring many-particle results to a phase space approximation. One observes 
that the number fluctuations AL^ are strongly increased during the sweep, 
but relax to a smaller value again afterwards. This evolution is well de- 
scribed within the semiclassical phase space picture. A further quantitative 
analysis of number squeezing during a Landau-Zener sweep is provided in fig- 
ure [71 comparing exact results (red) to an ensemble simulation (green). For 
a pure BEC with a given particle density, number fluctuations are given by 
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Figure 7: Number squeezing during a Landau-Zener sweep. Many-particle 
results (red) are compared to an ensemble simulation (green), (a) Evolution 
of the relative squeezing parameter £jy for a slow sweep with a = 0.1. (b) 
Final value of after the sweep as a function of the parameter velocity a. 
(c) Evolution of the spectroscopic squeezing parameter £| for a slow sweep 
with a = 0.1. (d) Dependence of the revival time of the occurrence of a 
squeezed state on the particle number N. The remaining parameters are 
J = 1, g = —5 and N = 200 particles. 



^^Ircf = (^2) / N . Thus one can define the parameter £jy = AL^/AL^ ref , 



which measures the suppression of number fluctuations in comparison to a 
pure BEC. Figure [7J (a) shows the value of £jy during a slow Landau-Zener 
sweep with a = 0.1. Indeed, £jy drops well below one for long times indi- 
cating number squeezing. Again, this feature is well reproduced by a phase 
space simulation (green). The final value of £jy after the sweep is shown in 
Figure [7] (b) as a function of the parameter velocity a. Number squeezing 
with £jy < 1 is observed for small values of a in the regime of the breakdown 
of adiabaticity, e.g. for large interaction strength, 2 J < g . The phase space 
simulation overestimates the variances and thus also £jy, but gives the cor- 
rect overall behaviour. For fast sweep, £jy tends to one as the state remains 
approximately coherent. 

However, an application of number squeezing in quantum metrology re- 
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quires a reduction of number fluctuations as well as a large phase coherence. 
Thus, a quantum state is defined to be spectroscopically squeezed if and only 
if 

A f 2 

ff : = N s < 1. (18) 

(i.> a + <iy> 2 

Spectroscopic squeezing indicates multipartite entanglement of the trapped 
atoms [301 HI]- The evolution of the squeezing parameter £| during a slow 
Landau-Zener sweep with a = 0.1 is plotted in figure [7J (c). While the 
number fluctuations AL 2 Z assume a small constant value after the sweep, the 
phase coherence (L x ) 2 + (L y ) 2 strongly oscillates due to the periodic de- and 
rephasing of the many-particle eigenstates (cf. figure 13]). True spin squeezing 
with £| < 1 is present only temporarily in the periods of maximum phase 
coherence. The timescale of the occurence of these minima depends linearly 
on the particle number N, as shown in figure [7] (d). For macroscopic particle 
numbers it takes very long before the states rephase such that £| < 1 is 
observed. Moreover these revivals are extremely sensitive to phase noise. 
Thus, it is doubtful that for realistic particle numbers Landau-Zener sweeps 
may be useful to generate squeezed states in a controlled way. Finally we 
note that the revivals of the phase coherence are not described by the phase 
space picture. Even small fluctuations in the phase coherence lead to large 
errors. Therefore the phase space approximation cannot account for the short 
periods where true spin squeezing £f < 1 is observed. 

Let us finally investigate the global dependence of the Landau-Zener tun- 
neling rate on the interaction strength g = UN in more detail. To this end 
we calculate the quantum and the classical tunneling rates given by equa- 
tions © and ffTOl . respectively, as well as the eigenvalues of the SPDM f|T5|) . 
We consider an initial state that is localized in the upper level for t — > — oo 
so that adiabaticity breaks down for a repulsive nonlinearity g > 2 J. As dis- 
cussed above, a change of the sign of the interaction strength g corresponds to 
an interchange of the two modes. For an attractive nonlinearity, adiabaticity 
breaks down in the lower level instead. Thus we obtain a global picture of 
the dynamics either by calculating the tunneling rate in the upper and the 
lower level for g > 0, or by calculating the tunneling rate in the upper level 
alone for g > and g < 0. In the following we choose the latter option. 

Figure [8] shows the results for J = 1 and a = 0.1, where the linear system 
evolves completely adiabatically. The left-hand side shows the many-particle 
and mean-field Landau-Zener tunneling probabilities as defined in equation 
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Figure 8: Left panel: Landau-Zener tunneling probability as a function of 
of the interaction strength UN for a slow parameter variation a = 0.1 and 
J = 1, N = 50. Right panel: Eigenvalues of the SPDM for t — > +oo. 
The exact many-particle results (red line) are compared to a phase space 
ensemble simulation (green diamonds) and the single-trajectory mean-field 
results (blue line). 

()6]) and (JTDjl . respectively. The right-hand side shows the eigenvalues of the 
SPDM ( TT5"]) for t — > +oo. Note, however, that the eigenvalues of the SPDM 
oscillate for t > as shown in figure SJ indicating a periodic loss and revival 
of coherence. Figure |S] shows the eigenvalues of the SPDM for large times, 
omitting the temporal revivals explicitly. As expected, adiabaticity breaks 
down as soon as g > 2 J and the Landau-Zener tunneling rate increases with 
g. In the adiabatic regime, one eigenvalue of the SPDM is close to zero, 
indicating a fully coherent state. Coherence is lost when the adiabaticity 
breaks down and particles are scattered out of the condensate mode. 

Figure [9] shows the results for a fast sweep a = 10 for J = 1. In the linear 
case, equation fll]) predicts a Landau-Zener tunneling rate of Plz = 0.7304. 
Surprisingly, the basic structure of the numerical results is very similar to 
the adiabatic case shown in figure |HJ The curves are shifted, but the general 
progression remains the same. This is understood as follows. As argued 
above, an attractive nonlinearity flattens the upper level so that Landau- 
Zener tunneling is decreased. The current example shows that this effect is 
so strong that the tunneling process is completely suppressed so that Plz — > 
for large negative values of g. On the contrary, a repulsive nonlinearity leads 
to an increase of Plz- The transition between an effectively adiabatic and 
non-adiabatic dynamics occurs at g = 2J for a slow parameter variation 
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a — > 0. For a fast sweep, Plz is non-zero in the linear case g = 0. However, a 
strong attractive nonlinearity can flatten the level so much that adiabaticity 
is restored again. Thus one can always enforce an adiabatic transition, but 
the necessary interaction strength \g\ increases monotonically with a. This 
behaviour is also reflected in the coherence properties of the final state shown 
on the right-hand side of figure [S] and EE 

One astonishing feature observed in the figures [S] and is the excellent 
agreement of the Landau-Zener tunneling rate Plz and the eigenvalues of 
the SPDM. Deviations are only found around g = in figure [9j This can be 
understood by a loss of the coherence between the two modes for long times, 
i.e. 



if we do not take into account for the temporal revivals illustrated in figure 
HI This happens either if the atoms are not in a coherent state any longer or 
if all atoms are localized in one of the modes. In any case we can rewrite the 
reduced SPDM as 



So the eigenvalues of the SPDM are directly given by the Landau-Zener 
tunneling rate if the two modes are not coherent. For strong nonlinear it ies g 
this is always the case and so the left- and the right-hand sides of the figures 
[8] and [9] show an excellent agreement except for a small region around g = 
in figure [9j In the non-interacting case (g = 0) the dynamics of all atoms 
is identical and the condensate will be fully coherent at all times. Thus the 
leading eigenvalue of the SPDM is always equal to unity independent of the 
Landau-Zener tunneling rate, such that the approximation f[2"0~j) is no longer 
valid in the non-interacting case. 

5 Semiclassical and adiabatic limit 

Having discussed various aspects of the mean-field many-particle correspon- 
dence in the previous sections, we now investigate the convergence to the 
mean-field limit quantitatively. The left-hand side of figure [10] compares 



(d{a 2 ) — > for i-> +oo, 



(19) 




(20) 
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Figure 9: Left panel: Landau-Zener tunneling probability as a function of 
of the interaction strength g = UN for a fast parameter sweep a = 10 and 
J = 1, N — 50. Right panel: Eigenvalues of the SPDM for t — > +oo. 
The exact many-particle results (red line) are compared to a phase space 
ensemble simulation (green diamonds) and the single-trajectory mean- field 
results (blue line). 

the mean- field Landau-Zener tunneling probability Plz(«) ( flOl) to the corre- 
sponding many-particle results (EJ) for different particle numbers and g = —5. 
While the many-particle dynamics usually converges rapidly to the mean- 
field limit, the occurence of a dynamical instability for \g\ > 2J leads to a 
breakdown of adiabaticity for small values of a. In this parameter region the 
convergence to the many-particle limit is logarithmically slow. This is further 
illustrated in figure [TD] on the right-hand side, where the Landau Zener tun- 
neling probability P^ p is plotted as a function of the inverse particle number 
1/JV. 

Another observation that can be drawn from the numerical data presented 
in figure [10] is that a simple mean-field description gives qualitatively wrong 
results in the adiabatic limit of small a. As already discussed in sectional the 
many-particle Landau-Zener tunneling probability P^(a) will always tend 
to zero for a — > since the level splittings in the many-particle spectrum may 
become small, but are always non-zero for finite N. Its mean-field counter- 
part (a) , however, is always affected by the appearance of the dynamical 
instability which destroys adiabaticity also for infinitesimally small values of 
a. Consequently, the Landau-Zener tunneling probability is believed to be 
non-zero even in this limit. This difference led to the claim that the adiabatic 
limit en — > and the semiclassical limit 1/N — > do not commute [31]. How- 
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Figure 10: Left panel: Landau-Zener tunneling probability Plz as a function 
of a for g = —5. Single-trajectory mean-field results (solid black line) are 
compared to exact many-particle results for N = 10 (green diamonds), N = 
20 (red squares) and N = 40 (blue circles) particles. The lines are drawn to 
guide the eye. Right panel: Slow convergence to the mean-field limit in the 
adiabatic regime. The Landau-Zener tunneling probability Plz is plotted as 
a function of the inverse particle number for a = 0.01 and g = —5. The 
black line corresponds to the single-trajectory mean-field results which are 
approached in the limit 1/N — > 0. 

ever, this claim is true only for the single-trajectory mean-field description 
which assumes a pure condensate at all times, which is obviously no longer 
true in the present case. 

As discussed in the previous section, the proper semiclassial limit of the 
quantum dynamics is a phase space flow rather than a single phase space 
trajectory. This description is valid also if the classical dynamics is unstable 
and the many-particle quantum state deviates from a pure condensate. The 
left-hand side of figure [11] shows the Landau-Zener tunneling probability 
P^ s (a) for different particle numbers calculated from the propagation of a 
semiclassical phase space ensemble as described in section HI It is observed 
that the many-particle results (cf. figure ITU]) can be reproduced to a very good 
approximation even for small values of a. Thus there is no incommutability 
of the adiabatic and semiclassical limits if the latter is interpreted correctly. 
Also the slow convergence to the single-trajectory limit is well described 
by the semiclassical phase space approach. The right-hand side of figure 
ITT1 shows -P L c z s as a function of the inverse particle number 1/N for a = 
0.01 which is well in the adiabatic regime. Significant differences to the 
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Figure 11: Left panel: Landau-Zener tunneling probability P{^ s as a function 
of a calculated from a semi-classical ensemble simulation for g = —5 and 
N = 10 (green diamonds), N = 20 (red squares) and N = 40 (blue circles) 
particles. The lines are drawn to guide the eye. The black line corresponds 
to the single-trajectory mean-field results which are approached in the limit 
1/N 0. Right panel: The Landau-Zener tunneling probability P{^ s as 
a function of the inverse particle number calculated from a semi-classical 
ensemble simulation for a = 0.01 and g = —5. 



many-particle results (cf. figure [TO]) are observed only for very small particle 
numbers, N < 10. 



6 Influence of phase noise 

We finally want to approach the question how an interaction with the en- 
vironment affects the transition from quantum-many body to the classical 
mean-field dynamics. To this end we consider the Landau-Zener problem 
subject to phase noise, which is the dominant influence of the environment 
provided that the two condensate modes are held in sufficiently deep trap- 
ping potentials [35l[36]. The many-particle dynamics is then given by the 
master equation 

= -i[H, p)~IJ2 ( h2 jp + - 2 ^i) • ( 21 ) 

U 3=1,2 

The effect of phase noise can be included in a single-trajectory mean-field 
limit starting from the dynamics of the Bloch vector (jSJ), whose evolution 
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equations are then given by [3"7H3"5] 



s 



s 



X 



y 



-2es y - 2Us y s z - 7s a 
2 Js z + 2es x + 2U s x s z 



(22) 



Thus, phase noise leads to transversal relaxation degrading the coherences s x 
and s y of the two condensate modes. Note that the magnitude of the Bloch 
vector is no longer conserved because of this effect. 

The resulting Landau-Zener tunneling probabilities are plotted in figure 
[T2las a function of a for different values of the noise strength 7. It is observed 
that phase noise has an important effect only for small values of a, where 
it drives the system to a completely mixed state with equal population in 
both wells so that Plz = 1/2. On the contrary, almost no consequences are 
observed for fast parameter sweeps. In this case, the tunneling time during 
which the atoms are delocalized is so short that phase noise cannot affect 
the dynamics. The transition to the incoherent regime occurs when the time 
scale of the noise 7 _1 is smaller than the tunneling time which is roughly 
given by a -1 . Therefore the sweep is incoherent such that Plz = 1/2 if 



while the interaction strength g has a minor effect only. 

Comparing mean-field and many-particle results, significant differences 
are observed for very small values of a and g = — 5 in the non-dissipative 
case 7 = 0, which has been discussed in detail above. In addition we note that 
already a small amount of phase noise is sufficient to remove these differences. 
For a — > and 7 / the mean-field approximation (122]) correctly predicts 
the transition to a completely mixed state with Plz = 1/2. Furthermore, 
significant differences are observed for g = —5 and intermediate values of 
a. In this case the many-particle quantum state is no longer a pure BEC 
but rather strongly number squeezed as discussed above. This state is more 
easily driven to a completely mixed state by phase noise than a pure BEC, a 
process which certainly cannnot be described by the simple single-trajectory 
mean-field approximation. 

Finally, these results suggest that Landau-Zener sweeps may actually 
be used as a probe of decoherence in systems of ultracold atoms (cf. also 
[3]). A measurement of the transition point to the incoherent regime where 
Plz = 1/2 gives an accurate quantitative estimate of the noise strength 7 
with a fairly simple experiment. 



oc < 7, 



(23) 
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Figure 12: Landau-Zener tunneling probability Pl,z{&) i n the presence of 
phase noise for g = — 1 (left) and g = —5 (right). The strength of the phase 
noise is chosen as 7 = (green diamonds), 7 = 0.01 (red squares) and 7 = 0.1 
(blue circles). Mean-field results (solid lines) are compared to many-particle 
results N = 40 particles (symbols). 

7 Conclusion and Outlook 

In the present paper we have presented an analysis of nonlinear Landau-Zener 
tunneling between two modes in quantum phase space. It was shown that 
adiabaticity breaks down if the interaction strength g = UN exceeds the crit- 
ical value 2 J- the Landau-Zener tunneling probability does not vanish even 
for an extremely slow variation of the system parameter. This phenomenon 
can be understood by the disappearance of adiabatic eigenstates in an inverse 
bifurcation in the mean-field approximation. Within the full many-particle 
description, the breakdown of adiabaticity results from the occurrence of di- 
abatic avoided crossings, where the level separation vanishes exponentially 
with the number of particles. 

The correspondence of the quantum dynamics and the 'classical' mean- 
field approximation has been discussed in detail. The many-particle and the 
mean- field Landau-Zener tunneling probability show an excellent agreement, 
because quantum fluctuations of the populations are small. In contrast, 
there is no fixed phase relation between the two modes, which certainly goes 
beyond the simple Bogoliubov mean-field theory. An improved classical ap- 
proximation using phase space ensembles can describe the depletion of the 
condensate mode and the loss of phase coherence as well as number squeez- 
ing £jy of the final state. Yet temporal revivals of this coherence are genuine 
many-particle effects and cannot be described classically. Thus, the spectro- 
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scopically relevant squeezing parameter £| is not reproduced by the ensemble 
simulation. However, the timescale for the occurence of these revivals and 
accordingly of the spectroscopical squeezing depends linearly on the particle 
number. For realistic setups, this is way to long compared to decoherence 
and phase noise rates. Before the sytem reaches the squeezed state, nearly 
all coherences are already lost. 

In the last section, we have studied how the dynamics depends on the 
number of particles N and compare our results to the discrete Gross-Pitaevskii 
equation that describes the dynamics in the limit iV — » oo. We show that 
the contradiction between the mean-field prediction and the exact many- 
particle transition rate in the adiabatic regime is no longer present in the 
phase space approach, and must therefore be considered as an artifact of 
the single-trajectory description. These results demonstrate the power of 
the phase space approach. However, in order to reproduce true quantum 
features such as quantum beats semiclassically, a more refined treatment is 
necessary. Semiclassical coherent state propagators have been studied in- 
tensively in single particle quantum mechanics in the limit h — > [12] • An 
extension to the mean-field limit of quantum many-body system must be 
based on the SU(M) phase space discussed in the present paper. In this 
case the particle number N is a number and not an operator and 1/N will 
serve as a proper semiclassical parameter. However, a numerical calculation 
for realistic particle numbers based on these methods, taking into account all 
relevant phase information between different trajectories, is as hard as the 
original quantum problem. 

Furthermore, we show that already the presence of a small amount of 
phase noise is sufficient to introduce enough decoherence to make the system 
'classical', so that the many-particle dynamics is well reproduced with a 
simple single-trajectory mean-field approach. Finally we have argued that a 
measurement of the transition to an incoherent Landau-Zener sweep could 
be used as a sensitive probe of decoherence. 
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